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New Theoretical Analysis of the LRRM Calibration 
Technique for Vector Network Analyzers 

Francesc Purroy and LIuis Pradell 



Abstract— In this paper, a new theoretical analysis of the four- 
standards line-reflect-reflect-match (LRRM) vector network-ana- 
lyzer (VNA) calibration technique is presented. As a result, it is 
shown that the reference-impedance (to which the LRRM calibra- 
tion is referred) cannot generally be defined whenever nonideal 
standards are used. Based on this consideration, a new algorithm 
to determine the on-wafer match standard is proposed that im- 
proves the LRRM calibration accuracy. Experimental verification 
of the new theory and algorithm using on-wafer calibrations up to 
40 GHz is given. 

Index Term— Calibration algorithms, calibration standards, 
calibration techniques, network-analyzer, on-wafer calibration, 
reference impedance, self-calibration. 



I. Introduction 

IT is well known that the systematic errors associated to a 
vector network-analyzer (VNA) can be modeled by error 
port adapters. The widely-used eight-term error model (Fig. 1) 
was first proposed in [1] and intended for two-way measure- 
ments. To compute its error terms, a calibration procedure is re- 
quired. A number of calibration techniques have been proposed 
in the literature. A technique that requires measurements from 
only three standards was first proposed in [2] (thru-short-delay, 
or TSD, technique). A more general approach (thru-reflect-line, 
or TRL, technique), in which the short standard is replaced by a 
highly reflective, but unknown, standard (reflect), was proposed 
in [3], and fully developed in [4]. Self-calibration techniques, 
which take full advantage of the redundancy in the calibration 
process, were developed in [5] and generalized in [6], [7] (Txx, 
Lxx). A particular case is line-reflect-match (LRM). Improve- 
ments in the computation algorithm of these techniques can be 
found in [8]. A variation of LRM (LRM-known Reflect) is de- 
scribed in [9]. 

The LRM technique requires two match (matched load) 
standards, one at each VNA port. As originally proposed in 
[ 6 L [7], both match standards should be equal and perfectly 
known. Since those requirements are not fulfilled in practice, 
the LRM calibration accuracy is reduced. The four-standards 
LRRM (line-reflect-reflect-match) calibration technique was 
proposed (but not mathematically developed) in [10], as an 
improvement with respect to LRM. The main advantage of 
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Fig. 1 . Error model for a vector network analyzer: (a) eight-term error model 
and (b) error matrices. 

o con. 



Fig. 2. Circuit model for on-wafer "match" standards. 

LRRM versus LRM is that only one match standard (partially 
unknown) is required to accurately compute the two unknown 
reflect standards T ri and r r2 , avoiding problems due to 
asymmetries between both match standards. The on-wafer 
match is modeled by a perfectly known resistor in series with 
an unknown inductance L e (see Fig. 2). The key point in the 
LRRM technique is an accurate determination of inductance L e 
from the computed reflection coefficient, r* lt of a measured 
reflect standard. In [10] a simple expression to compute L e 
from r* x is proposed [expr. (3) of [10]]. In this expression, 
it is implicitly assumed that a reference-impedance (to which 
the calibration is referred) does exist and equals the match 
impedance. However, this assumption is not justified, and its 
impact on the determination of L e is not considered. 

In this paper, a full theory for LRRM is developed. To the 
authors* knowledge, a rigorous analysis of LRRM has not been 
published yet. Based on this theory [14], it is shown that the 
'reference-impedance' (to which the LRRM calibration is re- 
ferred) cannot generally be defined whenever nonideal stan- 
dards are used. This consideration leads to a general expression 
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for r f c u that, in contrast with (3) of [10], does not assume a refer- 
ence-impedance. The expression is used to accurately determine 
L e , demonstrating accuracy improvements in the calibration. 

II. Theory of LRRM Self-Calibration Technique 

In this section, the three-standards self-calibration theory 
[6] is generalized to the four-standards case, in particular for 
the LRRM technique. It is assumed throughout this paper 
that the VNA systematic errors can be modeled by the bidi- 
rectional eight-term error model given in Fig. 1. The error 
terms are unknown and have to be determined through the 
calibration procedure. Four error terms (ebo»eoi»eio,eu) 
model the systematic errors corresponding to port 1 of the 
VNA (error-box X\ while e 22 , e 2 3, e 32 , and e 33 model 
the systematic errors corresponding to port 2 of the VNA 
(error-box Y). It is assumed that isolation terms (terms that 
account for direct transmission between port 1 and port 
2) are very small in practice and can be neglected (i.e., 
e 02 — ^20 = ei2 = e 2 i — e 0 3 — e 30 = e 13 = e 3 i = 0). 
Error-boxes are described by their transmission matrices ( 
J-matrices) [X] 9 \Y], respectively, defined as 



\Y] — 1 • [Fl = — • ( 632 ' 623 ' 
e 32 e 32 V 



e x o * coi - eoo • cu eoo 
-en 



1 J 



- e 22 <e 33 e 22 J ^ 

where [X), (F) are normalized versions of [X], [Y], respec- 
tively. 

A. Equation-Systems for the Error Matrices [X] and [Y] 

First, the calibration standards used in the LRRM technique 
are defined and their corresponding transmission-matrices 
(T-matrices) are given. 

1 ) LINE — Transmission line with a known length, propaga- 
tion constant and characteristic impedance 



(2) 



2) REFLECT1 — Dual one-port device made up by two 
identical, isolated loads with a reflection coefficient T rl . 
Its magnitude is high but unknown and its phase must be 
known within ±7r/2 radians 



(3) 



where SJJ (= 0) is the transmission 5-parameter of the 
standard reflect 1. 
3) REFLECT!— Dual one-port device made up by two 
identical, isolated loads with a reflection coefficient r r2 . 
Its magnitude is high but unknown and its phase must be 
known within ±7r/2 radians 



P3) = 



Cr2 
°2l 



V-r r2 i J 



(4) 



where SJ? (= 0) is the transmission 5-parameter of the 
standard reflect2. 



To obtain independent measurements, a nonideal open- 
circuit is used as reflect 1 and a nonideal short-circuit as 
reflect2 . 

4) MATCH— One-port device with an unknown reflection 
coefficient r m . Its magnitude should be very small. As 
reflectl, 2 and match are nontransmission devices, it is 
possible to define four virtual standards combining the 
match (port 1 or 2) with reflectl (port 2 or 1) or reflect2 
(port 2 or 1). If, for example, the match standard is con- 
nected to port 1, the following two standards are defined 
MATCH-REFLECT1: 



\Ta\ - - . ( ~^ ri • r m r m \ 
1 4j " s 2 t 1 I -r rl i ) 



(5) 



where SJft'' 1 (= 0) is the transmission S-parameter of 
the standard match-reflect 1. 
MATCH-REFLECT2: 



\ta — \ ( — ^1-2 • r m . r m \ 
lioJ ~s™* \ -r r2 i ) 



(6) 



where S£} r2 (= 0) is the transmission S-parameter of 
the standard match-reflect2. 
For each standard [2<J, a measurement transmission-matrix 
Mi (i ss 1, 2, 3, 4, 5) is defined as (see (A3) in Appendix A) 



M) = F)-[Ti\-F] 



1 



eioe 32 



(7) 



Combining measurement matrix [M x ] from standard 1 (line), 
with measurement matrix [Mi] from another standard i (i = 
2,3,4,5), the following linear-equation systems for matrices 
[X] and [F] are obtained from (7) 



[Mi] . [M x )~ l 
[Mi] • [Mi]" 1 



ppi-m-pii-py- 1 

(i =2,3,4,5) 



(8) 



or, in a more compact form 

[Mn].pn=pnp;i] 

[W«].lP)«[7).pli] 

(1=2,3,4,5) (9) 
where [Ma] and [T n ], are defined as 

[MiJ^MJ-iMJ- 1 (i = 2,3,4,5) (10) 
[Tn] —[Ti] • pi]"" 1 (i = 2,3,4,5). (li) 

In (8HH). matrices [Af<] (i = 1,2,3,4,5) are obtained 
from the measured (S)-parameters of standards [see (A5) in 
Appendix A], matrices (!•] (» = 1,2,3,4,5), partially un- 
known, are given by (2H6). and p?) and (F) are the unknown 
error-matrices. Equation (9) express a similarity transformation 
between matrices [Mu] and [Tn], with the two following 
mathematical properties. 

1) Trace conservation: 



2& l + Iff = + 



(12) 
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2) Determinant conservation: 



det([r«]) = det([M 41 ]). 



(13) 



Substituting (2)-{6) into (1 1) and the result into (12), we ob- 
tain 

1 _ M£+M% 
S r 2 \ (-e^.r^ + e-rU) 



u 2l - if! -re 

1 _ MH + 

1 _ M\\ + MH 

5 2 r l (-^• | »-r rl -r w + e-Ho 

1 MU + M?l 

(-er«..r p2 r jn + e-'r'0' 



(14) 



Equations (14) allow normalizing [Tu] and [Mn] to the same 
normalization factor to avoid the singularities of nontransmis- 
sion standards (S£ = = S£ rl = 5£ r2 = 0) in (3H7). In 
fact, using the normalized matrices [J54<] and [7<] (i = 2, 3, 4, 5) 
defined in Appendix A (A8), equation-systems (9) are written 
as 



pBal • (F) =(F) . psi 



where 



[Ma] =[5^1 • [Mx]" 

psi 



(15) 

(16) 
(17) 

(18) 
(19) 



ill =Pi|-pi] * 

Ti-_M±M_/'-r?i r rl \ 

2J (_rfH.I*+e-r') V-r rX 1 J 
Mg + Mjj (-T* 2 r r2 \ 

en- J*E±m% (-^nr m r m \ 
1 4J ~(-<^ • r rl • r m + e-7-o ^ -r rl 1 ; 

(20) 

— r^rm r?n \ 
1 oJ (- e yi-T r2 -r m + e-y') \ -r r2 1 

(21) 

Note that, since [Ma] in (16) do not require the measurement of 
^21* > the y are not singular; therefore, (Tu ] ( 1 7H2 1 ) are not sin- 
gular. Equation system (15) can be rearranged to form a system 
of 4 linear-equations for the 3 elements of pf ] 



(TX-MX _T?{_ -M$ _0_ \ 

T ii nt-MX 0 -Mg 

-Mj _0_ W-Ma _Jn__ 

0 -M$ t™-m??J 

/^iT\ /o 



"IS 
V i / 



(22) 



A similar system of 4 linear-equations for the 3 elements of 
[F] could also be written substituting X u , X^, X^* for 7Tu 
^L2, *2i in (22). Since 4 standards (i = 2, 3, 4, 5) can be com- 
bined with standard 1 (thru), equation-system (22) provides 16 
linear equations for [X] and [using TTT, W\ in (22)], 16 
linear equations for [F]. It can be shown (numerically) that the 
system composed by 16 equations has a maximum rank of 3. 
Therefore, only 3 terms of matrix [X] can be determined, from 
which error- terms e 0 o, en, ei 0 • eoi are computed. The same 
holds for [F] and e 22 , e 33 , e 2 3 • e 32 . The other two error- terms 
required (e l0 e 3 2) eoi^z) are computed from the measured ma- 
trix [Mi] (standard Thru') and the results obtained for e 00 , en, 
eio - coi and e 22 , e 3 3 t C23 • e 32 . 

There are 560 possible systems of 3 equations in (22), from 
which 16 have rank 2, The rest (544) have rank 3. The criterion 
proposed to select a three-equation system for the determina- 
tion of error-terms, is the system condition-number (CN). Using 
experimental calibration data, the CN for every system of rank 
3 in (22) has been computed. The experimental set-up is com- 
posed by a Hewlett-Packard 851 OB VNA and a Cascade-Mi- 
crotech SUMMIT 9000 wafer-probe station, with its calibration 
substrate type LRM-ISS. Measurement frequency-range was 1 
to 40 GHz. Two possibilities for connecting the match standard 
have been considered, match connected to calibration Port 2, 
and match connected to calibration Port 1. The following two 
best conditioned systems have been found. 
Match standard connected to Port 1 : 



-Ml 

-MU 



0 
0 



-Mil 

-nt 



m-Mii)[xr 2 \ 



MU -III. 



TJI_ 
-MU 

m 



IS - MU 



MJI 

- 1 -m 



m-Mii 



l 21 



-Ml? 



(I) 



Mji 

Match standard connected to Port 2: 



o -Wi 



(23) 



[ 51 

M$ 



51 




MH-Tg 
Ml 



-Ml 
-MU 



0 
0 

-Mil 




-7f 



MU- 



(24) 
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B. Expressions for T rl (Reflect 1) and T r2 (Reflect2) 

In this subsection, we focus on the determination of the un- 
known parameters T vX (reflectl) and T r2 (reflect2), in (3M<5), 
which are needed prior to solving systems (23) or (24). Using 
matrices [T 2 ] and [r 3 ] of standards reflectl and reflect2 [(3), (4), 
respectively], two new combined normalized matrices, [2^] = 
[U] ■ [T^]- 1 and [M32] = [M^\ • pEE)"" 1 , are defined. They can 
also be expressed in terms of p5i) f \^2i] and (3S?5i), pWji) as 



PS) ■ f^ 1 



(25) 
(26) 



Since these new matrices also fulfill trace and determinant 
conservation, (12) and (13) can be applied to (25) and (26) to 
obtain the following expression 



21 



T\2 , 7-> 2 l 



°31 



(27) 



where 

cll = Wi-Wx-Wx-Wi 

Substituting (1 8) and (19) into (27) and then using (12) again, 
the following equation for IVi and I\. 2 is obtained 



i% • (-i + Kit ■ (i - r* • e 2 ••>•<)) + r r2 . (2 • r rl ) 

, r>2 / 7x21 , \ r ^21 .- 2 .-vf 



+1*! • [Kl\ - 1) - J® • e- 



(28) 



where 



^21 _ 
A 3i — 



^31 



Repeating the above procedure for (T^) = pslj . [In]"" 1 and 
[243] = * p3i)~\ to following expressions are obtained 
for r r2 and r ri , respectively 



r r2 = 



where 



• (1 - rvi • r m • e*T> . Kit) + r m • (k$ - 1) 

(29) 



(j«+ji^-(Hff+Bff) 

r _ *1r • (g- 2n "' - r 2 2 ) + iVg . (r r2 - r m ) 
pl r,. 2 • (1 - r r2 • r m • 8 *n . k%) + r m • (^ - 1) 



(30) 



where 



^31 



(tf 4 \ l +*ff)-(w#+«2?) 



Substituting (29) [or (30)] into (28), a sixth degree equation 
for r ri (for T r2 ) is derived. It can be shown that it has the fol- 
lowing two double roots 



r rl = ±e-v* (r r2 = ± e -^) . 



(31) 



The other two roots are computed solving the following second- 
degree equation 



(r rl ) 2 -a + r rl .6 + c = o 



(32) 



where 



a • <T 2 -' • (2 • Kl\ - l) - • F 2 , ■ K\\ 

0 

21\ 2> 



r 2 



b =2 ■ I* • Kl\ . - 2 • JKg 4- j 

c . JKg . e- 2 ^' (2 . Kil - (Kit) 2 ~ l) 

Equation (32) is used to compute T rl . It can also be used to 
compute T r2 by exchanging r ri with T r2 , and Kl\ with K\\. 
To select the right root of (32) the phase of T rl (or r r2 ) must 
be known within ±tt/2 rad. Note that, prior to solving (32), 
T m should be known. An iterative method to compute T m is 
presented in the next subsection. 

C. Determination of T m (Match) 

According to [10], an on- wafer match standard can be mod- 
eled using the simple circuit of Fig. 2. Since the equivalent in- 
ductance L e is unknown, an initial estimation is used for T m , 
namely, T m = 0. This is equivalent to assuming L e = 0 
and R ~ Zq, where Zq is the normalizing impedance (usually 
50 ft). Other values could be used for R> provided they are well 
known. When using on-wafer standards, some match elements 
are trimmed to a known accuracy (typically better than ±2%). 
Using T m = 0, (32) is solved for an initial estimation of T rl 
(or T r2 ), r* x , T* 2 , respectively. The next step is to compute the 
actual value of r m by using T' rl (or r« 2 ). In [10], an expression 
is given that relates the computed (measured) reflect admittance 
(Y£) to its actual value (Yg) whenever the match standard is 
improperly defined (expression (1) of [10]) 



(33) 



where Y& » l/Z* = + and Yjg » 1/R are 

the actual and defined values for the match admittance, respec- 
tively. Equation (33) can be interpreted as a change in the cali- 
bration reference-impedance, from R (defined value) to Zfa = 
R -f jfwL 6 . In fact, substituting (33) into the expression for the 
computed reflection coefficient of standard reflect T* 



p 



(34) 
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where Y 0 = Yg 
obtains 



l/R is the normalizing admittance, one If t = 0, then (38) reduces to 



(35) 



Expression (-J5) can be interpreted as a change in the reference 
impedance to which an imperfect calibration is referred; 
imperfect means that some assumptions about the standards 
do not hold, in particular the actual and defined values for the 
match admittance are different. Therefore, T£ (measured by 
the imperfect calibration) is referred to the equivalent reference 
impedance^. The concept of reference impedance has been 
used in the literature to compare different calibrations [12], 
[13]. Obviously, a normalization impedance (Z 0 )> to which 
the actual (or assumed) reflection coefficients are referred, 
does exist. As we show later in this subsection, a reference 
impedance may not be defined for some imperfect calibrations. 
An expression equivalent to (35) is [1 1] 



r i-r M r* 



(36) 



where r„ - % + and T H = 
(Y 0 - Yg) I (Y 0 + Yg) . In other words, (33) assumes 
that a reference impedance does exist and it is equal to the 
match impedance, Zfa. To compute L e , in [10] it is assumed 
that the real part of Yg in (33) is zero (this assumption is 
reasonable for on-wafer open standards). In this case, the 
following expression [(3) of [10]] is obtained from (33) 



' (37) 



where G% and B% are real and imaginary parts of Yg 9 respec- 
tively. 

However, in an actual (nonideal) LRRM calibration the con- 
cept of calibration reference impedance may not be defined. 
Therefore, (33) and (37) cannot generally be applied. In fact, 
assume the following nonideal standards: 

1) Nonideal line with (nonzero) physical length (£ ^ 0) and 
a characteristic impedance Zl which is different from the 
normalizing impedance Z 0 . A reflection coefficient pl is 
defined for Z L , p L = (Z L - Z 0 )/{Z L + Z 0 ) t different 
from zero. However, the LRRM algorithm assumes a per- 
fectly matched line (pt = 0). 

2) Symmetrical unknown reflect (r rl ). 

3) Nonideal symmetrical match (T m ^ 0), but assumed 
ideal in the LRRM algorithm (T m = 0). 

As shown in Appendix B, the computed reflect reflection co- 
efficient with this actual LRRM calibration is given by (B.4) 
[see equation (38) at the bottom of the page]. 



(£m - r rl ) 
r rl -r m -i 



(39) 



Expression (38) shows that the actual reflection coefficient 
(T H ) and the estimated (T$ x ) values are not related by a change 
in the calibration reference impedance. Therefore, a reference 
impedance, to which the reflection coefficients computed by 
this imperfect calibration are referred, cannot be defined. In the 
particular case of zero-length line (i = 0), (39) holds, and the 
calibration reference impedance does exist and equals the actual 
match impedance Zfa, 

The equivalent inductance in the match model can be com- 
puted using (38) assuming that the line standard is perfectly 
matched (p L = 0). In this case, (38) reduces to 



• rl 



rl 



"e^-'.r m .r rl -i 



where 



_ (w<L e ) 2 + j-2-R.(wL e ) 
4-R 2 + (w<L e )* • 



(40) 



(41) 



Substituting (41) into (40) (using the negative sign) and 
enforcing |r rl | = 1 (reflect assumed fully reactive [10]), a 
second-order equation for (w • L e ) is obtained 

a • (w • L e ) 2 + 6 . (w . L e ) + c = 0. (42) 

where 

a =2 • »e(rj L ) + |r T %| 2 - 2 . SRe {T rl . e 2 ^') 

b =4 • R dc • (SJrn^J + S*m(rj x ■ c 2 ^ 1 )) 
c=4./i 2 c .(|r^| 2 -l). 

Equation (42) gives two solutions for the equivalent induc- 
tance I e . Equation (41) is used to select the solution giving the 
smallest |r m |. 

III. Simulations 

To show the advantages of this new method, it is useful to 
simulate the LRRM calibration, and to compare results to those 
of the method in [10]. Two cases are considered, ideal line, and 
nonideal line, 
a) Ideal Line 
Line: Perfectly matched, 1 ps delay, 
Reflectl: Symmetrical open-circuit, C - -12 fF. 
Reflect2: Symmetrical short-circuit, I = 6.244 pH. 
Match: R = 50 ft and L e - -7 pH. 
Using the proposed calibration algorithm, L e is com- 
puted and compared to the assumed value (L e = -7 pH), 
as shown in Fig. 3(a). Whereas the new method is giving 



r rl =± 



(r m -r,^.(o, L ) 2 -i) 



e 2 ^ ■ (pl - Tri) • (pl - r m ) + (r rl • *l - 1) . (i - r m ■ PL ) 



(38) 



1. 
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Fig. 3. Equivalent inductance (L c ) of a simulated "match" standard using the 
new method proposed in this paper (-) and the method proposed in [10] (-): 
(a) "line" perfectly matched and(b) "line" mismatched at low frequency (actual 
case). 



the right value, the method in [10] is giving L e ranging 
from -12 pH to -13 pH as a function of frequency. It can 
also be shown that the error increases when the assumed 
| X e | increases, 
b) Nonideal Line 
Line: Coplanar-waveguide transmission line, with a 
delay of 1 ps. Its characteristic impedance, assumed 
50 fi in the upper frequency range, was modeled using 
the expression for a lossy line with the following 
parameters 



7 a_ [R + 7 



w • L 



-w-C 

L =386 nH/m C=154pF/m 



G ~w • C - tg(5 L R oc 



where tg(5i,) is the dielectric (alumina) loss-tangent 
and <xq is the conductor (gold) conductivity. The scale 
factor for R was found by fitting (as a function of 
frequency) the loss-constant of the coplanar-waveguide 
line, that was obtained from a TRL calibration per- 
formed on-wafer. 

Reflectl: Symmetrical open-circuit, C = -7.19 fF. 

Reflect2: Symmetrical short-circuit, L = 6.24 pH. 

Match: R = 50 Q and L e = -7 pH. 
The result of this second simulation is shown in Fig, 3(b). At 
low frequencies, where the line standard is not well-matched, 
both methods get bad results, because the algorithm is assuming 
that the line is perfectly matched. At higher frequencies the re- 
sults are similar to case (a). 

IV. Experimental Results 

The calibration algorithm proposed in this paper has been ex- 
perimentally tested using the experimental set-up described in 
the preceding section. Details of the calibration standards and 
wafer-probes are 

CALIBRATION SUBSTRATE: LRM-ISS (Cascade-Mi- 

crotech). 

Line: 1 ps delay line assumed perfectly matched. 
Reflectl: Open circuit (probe tips in air) assumed sym- 
metrical. 
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Fig. 4. Equivalent inductance (L c ) of a measured "match" standard computed 
using the new method (-) and the method in [ 1 0] (-). 
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Fig. 5. Measured |S U I of a 40 ps-delay open-ended coplanar-waveguide line 
using LRRM (-) and LRM (-). 



Reflect2: Short circuit assumed symmetrical. 

Match: Coplanar load, R = 50 fi. 
WAFER-PROBES: WPH-305-150 (Cascade-Microtech). 
Fig. 4 shows the equivalent inductance L e computed using the 
method proposed in [10] as well as using the new method pro- 
posed in this paper. Above 10 GHz, L e is more frequency- in- 
dependent when using the new method. In the low-frequency 
margin, both methods compute incorrect values, because the 
line standard is not well matched. These results are in agree- 
ment with predictions of simulations described in Section III 
(see Fig. 3). 

Fig. 5 shows the reflection coefficient magnitude of a 
40 ps-delay open-ended coplanar-waveguide line (not used 
for calibration), measured using LRM as well as LRRM. It 
can be seen that results are very similar. Therefore, the new 
LRRM theory and calibration algorithm are validated. The 
small differences observed are due to asymmetry in the match 
standards using LRM, Differences are more noticeable at 
high frequencies due to increasing asymmetry. Since LRRM 
requires only the measurement of a match in one of the two 
ports, significant accuracy improvements of LRRM calibration 
using the new algorithm, compared to LRM, are expected at 
higher frequencies (beyond 40 GHz). 

V. Conclusion 

In this paper a new theory of the LRRM calibration tech- 
nique, that generalizes the three-standard self-calibration theory 
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for VNA to four standards, has been presented. Using the crite- 
rion of equation-system condition-number, the optimum equa- 
tion-system has been selected from experimental on- wafer data, 
to compute the coefficients of the eight-terms VNA error model. 

A theoretical study on the reference impedance associated 
to an imperfect (actual) LRRM calibration has been developed, 
showing that the reference impedance is not defined whenever 
the line length is not zero. Therefore, a new method to com- 
pute the equivalent inductance L e of the match standard model 
is proposed. Simulations show that the new method accurately 
computes the assumed L e , in contrast with a former method pro- 
posed in the literature. 

Experimental on-wafer results demonstrate the feasibility 
of the new theoretical LRRM formulation. The measurement 
of a 40 ps delay open-ended line shows the advantages of the 
LRRM algorithm versus LRM whenever asymmetries in the 
match occur, in particular at high frequencies. 

Appendix I 

Definition of Self-Calibration Measurement Matrices 

Referring to Fig. 1, we have for the measured standard i (i = 
1,2,3,4,5) 

• In forward measurements: 

(S)"^" (J) (A,) 

• In reverse measurements: 

(5)-S5S«»>-™-""-tt) (A2) 

where [XJ, [F] are defined in (1). 

Combining (Al) and (A2), the measurement transmission 
matrix [Mi] is defined [6] 



(A3) 

Developing (A3), the 5-parameters corresponding to [AfJ, are 
expressed as 



°0' a 3 



1 " 



5a., 



1 - 



VQ3 

3'QO 



o 2 2 







1- 






a 3 -a 0 


I" 


a^fc 3 






1- 





with 



1 / 9 M 



10 *-'21 



CMi qMi qM{ \ 

fll °22 °11 \ 

i J 



(A4) 



(A5) 



A, B, C, D, and two additional measurements ( B, F), defined 
as: 



a 0 



a$ f ao 

B D -— F=4 
a 0 o ; 3 a' 3 



(A6) 



In case of dual one-port devices (reflect 1, reflect2, match- 
reflect 1, match-reflect2), a' 0 ~ b' 0 - a 3 = 63 = 0, and (A4) 
reduces to 



Normalized measurement matrices can also be defined from 
(A5) as [Ml] = [Mi] • S%[ 1 . In a similar way, normalized ma- 
trices of standards are defined as pi] « [T<] • 5^. Using 
the 5-parameter definition (A.4), the ANA internal switch is 
included in the measurement matrix and the eight-term error 
model holds for both directions. The number of required mea- 
surements to compute [Mi] from (A4) is six, four standard ratios 



Q>Q 

a 3 



(A7) 



and their normalized measurement matrices are (i = 2, 3, 4, 5) 

-si> 22 I 1 )- (A8) 



Appendix II 

Expression for the Computed Reflection Coefficient 
of Reflect Standard 

In this appendix, the relationship between the computed re- 
flection coefficient T* x of standard reflect and its actual value 
T rl is derived. The calibration algorithm assumes that standards 
match and line are perfectly matched, but they actually have a 
reflection coefficient different from zero (T m ^ 0 and p L £ 0) 
referred to a given normalization impedance Z 0 . Therefore, the 
actual line standard transmission matrix is 



pi] - J 



1 



(PL) 3 

e- 1 "' - (pl) 2 • e T ' 



p L ■ (e^ - ( 



c-Hj 



(Bl) 



Matrices [I2] to P5) [(18H21)] do not change. The normal- 
ized measurement matrices that VNA would measure can be 
computed using (IS) 



[M 21 ]=pT]-p ; 2i]m- 1 

p^iHm-pnaii-m- 1 



(B.2) 



where matrices P21], [?3i], [T 5 i] are computed using (17) (i = 
2, 3, 5) with pi] given by (B.l). Since the line and the match 
standards are assumed perfectly matched, the unknown reflec- 
tion coefficient of standard reflect is computed using (32) 
with r M = 0 



(B3) 



where the definitions of iff J , iff J are given in expressions (28) 
and (29). Substituting (28) and (29) into (B3), the following ex- 
pression is obtained [see (B4) at the top of the following page]. 



